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ABSTRACT 

We present nine near-infrared (NIR) spectra of supernova (SN) 2005cf at epochs from 
— 10 d to +42 d with respect to £?-band maximum, complementing the existing excel- 
lent data sets available for this prototypical Type la SN at other wavelengths. The 
spectra show a time evolution and spectral features characteristic of normal Type 
la SNe, as illustrated by a comparison with SNe 1999ee, 2002bo and 2003du. The 
broad-band spectral energy distribution (SED) of SN 2005cf is studied in combined 
ultraviolet (UV), optical and NIR spectra at five epochs between ~8d before and 
~10d after maximum light. We also present synthetic spectra of the hydrodynamic 
explosion model W7, which reproduce the key properties of SN 2005cf not only at 
UV-optical as previously reported, but also at NIR wavelengths. From the radiative- 
transfer calculations we infer that fluorescence is the driving mechanism that shapes 
the SED of SNe la. In particular, the NIR part of the spectrum is almost devoid of 
absorption features, and instead dominated by fluorescent emission of both iron-group 
material and intermediate- mass elements at pre-maximum epochs, and pure iron-group 
material after maximum light. A single P-Cygni feature of Mgll at early epochs and a 
series of relatively unblended Co II lines at late phases allow us to constrain the regions 
of the eject a in which the respective elements are abundant. 

Key words: supernovae: individual: SN 2005cf- supernovae: individual: SN 1999ee - 
supernovae: individual: SN 2002bo - supernovae: individual: SN 2003du - techniques: 
spectroscopic - radiative transfer 



1 INTRODUCTION 

An empirical relation betwe en their maxi mum brightness 
and light curve decline rate makes Type la 

supernovae (SNe la) one of the best to ols to measu r e cos - 
mol ogical distances. Using th is approach Rie ss et all ([l998) 
and IPerlmutter et al.l ([19991 ) discovered that our Universe 
undergoes an accelerated expansion. This was in turn in- 
terpr eted as evidence for the existence of 'dark energy' 
(e.g. iRiess et "all l2007h . Despite this importance for astro- 
physics and cosmology our knowledge about the progeni- 
tor systems or the explosion mechanisms of these luminous 
events is still limited. They are commonly believed to be 



thermonuclear explosions of carbon-oxygen white dwarfs in 
binary systems. However, many aspects are still unclear; 
e.g. the masses of the exploding stars, the nature of the 
companion star, the ignition conditions, the exact mech- 
anism of the burning process or the origin of the various 
subclasses of SNe la. Systematic uncertainties arising from 
this ignorance are one factor limiting the precision of cosmo- 
logical measurements. One way to improve our understand- 
ing of these objects is the comparison of observed spectra 
with predictions from theoretical models. This, however, re- 
quires complete da ta sets with a wide wavelength coverage 
(jRooke et al.ll2012l ). 

In particular, the near- infrarecfl (NIR) regime has been 



* Based on observations collected at ESO, Paranal. Program ID 
075.D-0823(B). 

f E-mail: egall@mpa-garching.mpg.de 



1 In this paper we will refer to wavelengths between 9500 and 
26 000 A when talking about 'near-infrared'. 
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proposed as a useful t ool to discriminate between different 
explosion models (e.g. M arion et al.ll2003h . This is mainly 
due to the fact that the lower opacity in the NIR enables 
us to look deeper into the eject a and thus probe layers still 
hidden to optical observations at the same epochs. Also, the 
NIR is believed to be less severely affected by line blend- 
ing compared to optical or ultrav i olet (UV) wavelengths 
(lHarknesslll99ll: IHoflich et al.lll993l : ISpvromilio et al.lll994l : 
IWheeler et al.l 19981 ). At the same time, however, observa- 
tions suggest SNe la to be exc e llent NIR standard can dles 
(lElias et al.l Il985l: iMeikld l2QQ0l : iKrisciunas et all I2QQ4I . but 
see lKattner et al.ll2012h , in contrast to optical wavelengths 
where they can only be used as standardzzea6/e candles. 
Thi s result was confirmed by a series of toy explosion mod- 
els ([Kasen 2006). While the homogeneity in the NIR might 
point towards a common explosion mechanism for all SNe la, 
it could also mean that the NIR photons are not as sensitive 
to variations in the eject a as previously suggested. 

While optical spectra have been studied quite 
extensively, the NIR regime is still comparatively un- 
explored. Although there are quite a few publications 



Table 1. Properties of SN 2005cf, as derived in the literature. 



19871: Sovromilio et al.1 I1992L 20041: Meikle et al.1 Il996l: 


Bowers et al.1 Il997l: lHamuv et al.1 2002allbl: IHoflich et al.1 


20021: iRudv et al.1 12002 


: iMarion et al.l 200l 120091: 


Benetti et al. 120041: lElias-] 


losa et al.l 120061; IStanishev et al.l 


20071: iTaubenbereer et al.l 


2011) one has to realize that for 



many SNe there is only one spectrum available - typically 
at post-maximum epochs - or that the objects are of a 
peculiar subclass. The number of well-sampled normal 
SNe la with pre- maxi mum- light ob s ervati ons in the NIR 
is actually quite small. iMeikle et al. Jl996l) have published 



six spectra of SN 199 4D. lHamuv et al.1 ( 1 2002 al lri) eleven 
spectra of SN 1999ee, iBenetti et al.l (|2004r ) six spectra of 
SN 2002bo. IPignata et al.l (l2QQ8h four spectra of SN 2002dj, 
Elias- Rosa et al.1 (l2QQ6h ten spectra of SN 2003cg and 
Marion et al.l (|2QQ9h four spectra of SN 2005am. Moreover, 
there are 14 spectra of SN 20 3du out to the nebular phase , 
pub lished bvlMotohara et al.l (l2006h , IStanishev et al.l (l2007h 
and lMarion et al.l (|2009h . 



In order to increase the sample of SNe la with good NIR 
coverage we present and analyze a set of nine NIR spectra 
of SN 2005cf, ranging from -10.2 d to +41.5 d with respect 
to_B-band maximum. SN 2005cf, according to I Wang et al.l 
(120091 ) ' the golden standard Type la supernova', is a spec- 
troscopically normal SN la with a more complete data set 
than most other SNe la. Optical and NIR photometry as 
well as o ptical spectroscopy o f SN 2005cf have been pub- 
lished bvlGaravini et all (j2007h . IPastorello et al.l (l2007h and 
IWang et al.l (l2009h . UV photomet r y and spectr oscopy have 
been published by iBufano et al.l ([2009) and I Wang et al.l 
(|2Q12h . 

Table [1] summarizes the most important properties of 
SN 2005cf: at ff- band maximum, it h ad a magnitude of 
A^max = -1 9.39 (IPastorello et al.ll2007h . AmisCB) lies be- 
tween 1.07 (IWang et al.1 l2009h and 1.12 (IPastorello et al.1 
l2007h. and the der ived 56 Ni mass is about .7 M^ 
(IPastorello et alJl2007h or 0.78 M (IWang et al.ll2009h . This 
corresponds quite well to the average values obtained for 
normal SNe la. An important shortcoming of the data set of 
SN 2005cf has so far been the lack of published NIR spec- 





Pastorello et al. 


Wang et al. 




(^2007) 


(2009^1 


JD Of Bmax 


2 453 534.0 ± 0.3 


2 453 533.66 ± 0.28 


ms max 


13.54 ± 0.02 


13.63 ± 0.02 


M s max 


-19.39 ± 0.33 




Ami 5 (B)true 


1.12 ± 0.03 


1.07 ± 0.03 


M( 56 Ni) 


O.7M 


(0.78 ± 0.10) M 


/i(MCG-01-39-003) 


(32.51 ± 0.33) mag 


(32.31 ± 0.11) mag 


e(b - v) host 


mag 


(0.10 ± 0.03) mag 


Table 2. Epochs of the NIR spectra of SN 2005cf. 


UT date 


Julian Day Epoch 


Instrument 


2005/06/02 


2 453 523.8 -10.2 d 


ISAAC 


2005/06/03 


2453 524.5 -9.5d 


NICS 


2005/06/07 


2453 528.6 -5.4d 


ISAAC 


2005/06/12 


2453 533.7 -0.3d 


ISAAC 


2005/06/14 


2453 535.8 +1.8d 


ISAAC 


2005/06/23 


2453 544.7 +10.7d 


ISAAC 


2005/06/26 


2 453 548.4 +14.4 d 


NICS 


2005/07/13 


2453 565.4 +31.4d 


NICS 


2005/07/24 


2 453 575.5 +41.5 d 


NICS 



ISAAC = 8.2 m Very Large Telescope + ISAAC, 

short- wavelength low-resolution spectroscopy (SWS1-LR); 

http: //www. eso.org/ sci / facilities / par anal / instruments / isaac / 

NICS = 3.58 m Telescopio Nazionale Galileo + NICS, Amici 
prism; http: // www.tng.iac.es / instruments /nics / 

troscopy. This deficiency is now remedied with the data pre- 
sented in this work. 

The paper is organised as follows. In Section [2] the NIR 
spectra of SN 2005cf are presented and the process of data 
reduction is described. In Section [3] we compare these ob- 
servations to the NIR spectra of other SNe la and construct 
combined UV-optical-NIR spectra to investigate the broad- 
band spectral energy distribution (SED). Sy nthetic spectra 
of th e hydrodynamic explosion model W7 ([Nomoto et al.l 
Il984h are used to explain the SED in terms of fluorescence 
and to perform a NIR line identification (Section 2}. Con- 
clusions are drawn in Section [5] 



2 OBSERVATIONS AND DATA REDUCTION 

The spe ctroscopic data of SN 2005cf were obtained with 
ISAAC dMoorwood et al.| [l999). mounted at the Very Large 
Telescope at the Paranal Observatory of ESO (European 
Southern Observatory), and NICS (Near-Infrared Camera 
Spectrometer) , mounted at the Telescopio Nazionale Galileo 
of Fundacion Galileo Galilei and INAF (Istituto Nazionale 
di Astrofisica) . The spectra were collected by the European 
Research Training Network (RTN) 'The Physics of Type la 
Supernova ExplosionsQ- They cover the time interval from 
10 d before to 42 d after 5-band maximum (see Table [2] for 
details) . 



2 http: //www. mpa-garching. mpg.de / ~::tn / 
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The spectra taken with ISAAC (at epochs —10. 2d, 
— 5.4 d, —0.3 d, 1.8 d and 10. 7 d with respect to S-band max- 
imum) were taken in the short-wavelength low-resolution 
setup (resolution R = & 500). The observational data 
are divided into the four bands SZ, J, SH and SK, which 
cover the wavelength range from 9000 A to 25 600 A. The 
NICS spectra (at epochs -9.5 d, 14.4 d, 31.4 d and 41.5 d 
with respect to S-band maximum) were taken with an Amici 
prism as disperser. They cover the whole NIR region between 
about 7500 A and 26 000 A at very low resolution (R & 50). 

For the NICS data, cross-talk and distortion corrections 
were applied with SNAF0. All other reduction steps for both 
the NICS and the ISAAC data were performed using irafQ, 
following standard procedures in the NIR. This included 
dark-current and flat-field corrections, and a pairwise sub- 
traction of sub-exposures taken with the target off-set along 
the slit in order to remove the sky emission. One-dimensional 
spectra were obtained from the two-dimensional ima ges us- 
ing a n optimal, variance- weighted extraction method (|Hornd 
1 19861 ) . For the ISAAC spectra the wavelength was calibrated 
with the help of arc-lamp spectra taken along with the SN 
exposures using exactly the same telescope configuration, 
and the calibration was checked with the position of known 
telluric features. For the NICS spectra, the wavelength cal- 
ibration was accomplished with tabulated values. Telluric 
features were removed from the SN spectra by dividing 
through spectra of telluric standard stars of type AO or simi- 
lar, which had been observed at a similar airmass as the SN. 
A first relative flux calibration was provided by multiplica- 
tion with a Vega spectrum. Then the absolute fluxes were 
adj usted to agree w i th co ntem poraneous JHK ph otometry 
bv IPastorello et al.l <j2007h and IWang etafl <|20Q9h . Finally, 
redshi ft (z 0.00646) and reddening \E(B - V) Ga i = 0.097 
mag, ([SchlegeL Finkbeiner. fe D avis 1998)] corrections were 
applied. 



3 RESULTS 

The nine final, reduced NIR spectra of SN 2005cf are pre- 
sented in Fig. [1] Before maximum light, large parts of the 
spectra resemble the infrared tail of a black-body contin- 
uum. This changes significantly at later epochs, when the 
spectra start to show increasingly pronounced spectral fea- 
tures. Around 10 d after maximum (and thereafter) two dis- 
tinct humps form between 15 000 A and 18 00 A. These have 
been previously identified as Co n emission ([ Wheeler et al.l 
1998). Characteristic features in the K band, which form at 
even later epochs, ha ve also been attributed to Co emission 
(Wh eeler et al.|[l998h . Compared to a corresponding black- 
body continuum there is an increasing flux deficit in the 
J-band region after maximum light. 

In the following we will take a more detailed look at the 
spectral properties of SN 2005cf, comparing the NIR data of 



3 Speedy Near-IR data Automatic reduction Pipeline, 
http: //www. arcetri. astro. it/ ~:ilippo/snap/ 

4 IRAF (Image Reduction and Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Figure 1. Time sequence of our NIR spectra of SN 2005cf. The 
gaps between the NIR bands are regions where the Earth's atmo- 
sphere is opaque. 



SN 2005cf to those of other well-observed normal SNe la and 
studying the SED in combined UV-through-NIR spectra. 

3.1 Comparison of SN 2005cf to other Type la 
supernovae 

Representative, normal SNe la, for which multi-epoch 
NIR spectroscopy is available, are e.g. SNe 1999ee 



(lHamuv et al 2002a, b) , 2002b o (jBenetti et al.l l2004h and 
2003du fStanishe v et al.l l2007h . An overview, comparing 
some of their properties to those of SN 2005cf, is given in 
Table El 

The light-curve decline rate Amis (B) of SN 1999ee 
characterizes it as a relatively slowly declining SN la. 
However, its peak magnitude agrees with t he light-curve 
shape as described by the Phillips relation ([Phillips! Il993l ; 
IPhillips et al.l Il999f). Marginally assign ed to the 'shallow- 
silicon' class bv branch et all (|2006h , the spectral fea- 
tures show that, ov erall, SN 1999ee is a typical SN la 
|Hamuv et al 1 l2002al lbh. Concerning the spectral features, 
the same can be said about SN 2002bo, though at early 
phases this SN has higher ejecta velocities than the aver- 
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Table 3. Comparison of characteristic properties of SNe 1999ee, 2002bo, 2003du and 2005cf. 





SN 1999ee 


SN 2002bo 


SN 2003du 


SN 2005cf 


Branch subtype" 


shallow-silicon 


broad-line 


core- normal 


core- normal 


Benetti subtype 6 


LVG 


HVG 


LVG 


LVG 


v (kms-M- 1 ) 6 


42 ± 5 


110 ± 7 


31 ± 5 


35 ± 5 


vio(Sin) (lO^ms" 1 ) 6 


9.90 ± 0.10 


11.73 ± 0.15 


10.10 ± 0.10 


9.94 ± 0.30 


M Binax 


-19.85 ±0.28 


-19.41 ±0.42 


-19.34 ±0.16 


-19.39 ±0.33 


Ami 5 (5) 


0.94 ± 0.06 


1.13 ± 0.05 


1.02 ± 0.05 


1.12 ± 0.03 


E(B - V) 


0.32 mag 


(0.43 ± 0.10) mag 


(0.01 ± 0.05) mag 


0.097 mag 


Redshift z 


0.0117 


0.0043 


0.0064 


0.00646 


References 


Hamuv et al. f2002a.b) 


Benetti et al. (2004, 2005) 


Stanishev et al. (2007) 


Garavini et al. (2007) 




Stritzinger et al. (2002) 


Stehle et al. (2005) 


Benetti et al. (2005) 


Pastorello et al. (2007) 


Benetti et al. (2005) 


Branch et al. (2006) 


Branch et al. (2009) 


Branch et al. (2009) 


Branch et al. (2006) 






Wane et al. (2009) 



a Sub-classification according to iBranch et al ] (120061 ). 

b Low-velocity-gradient (LVG) and high-velocity-gradient (HVG) SNe as defined by Ben etti et al ] (120051 ). 



age normal SN la (jBenetti et al.l 120041 ; IStehle et all l2005h 
and has been included in the 1 bro ad- line' and ' h i gh- vel ocity- 
grad ient' SN la su bclasses by IBranch et all (|2006l ) and 
IBenetti et all (120051 ). respectively. SN 2003du, finally, is very 
similar to SN 2005cf in many respects , both being 'core- 
normal', 'low- veloc ity- gradient' SNe la (Branch et al ] l2006l : 
IBenetti et al.ll2005l ). Fig. [2] shows that the NIR spectral evo- 
lution of all four SNe is remarkably similar, with a few ex- 
ceptions worth mentioning. 



3.1.1 Spectra ten to five days before B-band maximum 

Panel (a) of Fig. [2] shows spectra between 10.2 and 5.4 d 
before B-band maximum. 

Around 10 400 A there is a very prominent absorption 
feature in SN 200 2bo. This fea t ure w as first observed in 
SN 1994D, where iMeikle et all (jl996h suggested that it 
comes from ei t her M gn A10926 or Hei A10830. However, 
I Wheeler et al.l (| 19981 ) sh owed evidence that t he line should 
come entirely from Mg. IBenetti et al.l (|20 04) followed this 
reasoning and attributed the feature in SN 2002bo to Mgn. 
The same feature, though much weaker, is also present in the 
spectra of SNe 2003du and 2005cf, but only marginally de- 
tectable in SN 1999ee. Whether this variation is a sign of dif- 
ferent Mg abundances or just different ionization and excita- 
tion conditions is un clear. A feature aroun d 1 6 800 A can be 
seen in all four SNe. Ham uv et al ] (l2002allbl). Marion et al. 
(120031 ). IBenetti et all <l2004h and IStanishey et al I (|2007l ) 
agree that it is likely due to Si II. IBenetti et al. I <|2004h believe 
that a wea k feature in SN 2 002bo at about 20 800 A comes 
from Sin. IStanishev et all (2007) instead propose Si III as 
the origin of the same feature in SN 2003du. In SN 2005cf 
the feature can also be discerned in the — 10 d spectrum, but 
appears more prominent in the — 5d spectrum. 



3.1.2 Spectra around ten days after B-band maximum 

Panel (b) of Fig. [2] shows spectra between 8.5 and 12.4 d 
after B-band maximum. 

In contrast to the spectra before maximum light, the Mg 
and Si features are no longer present in the spectra 3 weeks 
later. Instead, the NIR spectra now show markedly more 



structure and are dominated by singly ionized Co and Fe. 
In particular, there are strong Co II features around 15 500 A 
and 17 500 A in all four SNe. 



3.1.3 Spectra around thirty days after B-band maximum 

Panel (c) of Fig. [2] shows spectra between 29.4 and 31.5 d 
after B-band maximum. 

The strong features in the H band seem almost iden- 
tical in the four Type la SNe, but especially similar in 
SNe 2002bo and 2005cf. In SN 2003du this region appears 
to be somewhat better resolved into individual, narrower 
features, which could point to a slightly lower maximum ex- 
pansion velocity of the iron-group-element-rich core of the 
eject a. Above 20 000 A, a characteristic pattern of Coil lines 
(Wh eeler et al. 19981) is now very prominent in all four SNe. 
As shown by iMarion et al.l (|2009h , these features appear in 
all SNe la at sufficiently advanced epochs. Similarly, a lack 
of flux in the J band, compared to a corresponding black- 
body continuum, is now clearly visible in all the SNe of our 
sample. 



3.1.4 Mg line velocities 

In previous work it was reported that the blueshift of the 
Mgn A10926 a bsorption linq j does not evolve with time 



lamu v et al 



Meikle et al 



in SN e 1994D (IMeikle et al.l 1 19961). 1999ee { 
l2002al lbh and 2003du (jStanishev et al.l l2007h 
([ 19961 ) interpreted this as the feature being detached, which 
means that it forms well above the region of peak emissivity 
at these wavelengths. Accordingly, the measured Mg veloc- 
ity will not change when the emission region recedes deeper 
into the core due to the increasing dilution of the eject a with 
time. Thus, the measured Mg 11 velocity traces a zone with 
high Mg 11 abundance, or specifically the lower boundary of 
such a zone where the density is highest. Mgn lines in the 



'Mgll A10926' is a short notation for the triplet 



A10914: 2p b 3d 2 D 
A10915: 2p 6 3d 2 D 
A10952: 2p 6 3d 2 D 



5/2 
3/2 
3/2 



2p b 4p 2 P 
2p 6 4p 2 P 
2p 6 4p 2 P 



3/2^ 
3/25 
1/2- 
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10000 15000 20000 25000 30000 
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10000 15000 20000 25000 30000 



Figure 2. SN 2005cf in comparison with SNe 1999ee, 2002bo and 2003du. Panel (a): about 12 to 5 d before B-band maximum. Panel 
(b): 9 to 12 d after S-band maximum. Panel (c): 29 to 32 d after £>-band maximum. 



optical part of the spectrum are less suited to unveil the 
location of the Mg-rich zone since i) they are intrinsically 
stronger so that very little Mg (already on the level of Mg 
present in the proge nitor star) is sufficient for them to form 
(jMarion et al.ll2003h , ii) they are more strongly blended with 
lines from other species, and iii) the emissivity at shorter 
wavelengths peaks at much larger radii than in the NIR (see 
Section T4. 2 p - at early phases inside or even above the zone 
of highest Mg n abundance. 

In order to verify whether the constancy of the Mgll 
A 10926 velocity may indeed be a generic feature of normal 
SNe la, we measured the position of this line in our ISAAC 
spectra of SN 2005cf. We obtain velocities of - 11 400 km s -1 
for the -10.2 d, -0.3 d and +1.8 d spectra, and -11200 
kins" 1 for the +5.4 d spectrum. Within the typical mea- 
surement uncertainties of —100-200 kins" 1 this is a per- 
fectly constant trend from the earliest spectrum to the dis- 
ap pearance of the line around maximum light, supporting 
the lMeikle et alJ (l996) interpretation of a detached feature. 
The absolute value of the Mgn A10926 velocity of SN 2005cf 
is also q uite similar to those derived for SNe 1994D (—11 700 
kms" 1 . iMeikle et al.1 [l99(jh and 1999ee (-10 500 kms" 1 , 
lHamuv et al.l l2002a,b), indicating a similar layering in these 
three low-velocity-gradient SNe. 

Unfortunately there is only a single pre-maximum NIR 
spectrum available for SN 2002bo, and accordingly the time 
evolution of the Mgll A10926 velocity cannot be checked. 
The value of — 14 400 kms -1 measured in the — 8.5d spec- 
trum, however, is significantly larger than in the other SNe. 
SN 2002bo is a high-velocity-gradient SN with particularly 
broad and strongly blueshifted Sin and Call lines in the 
early optical spectra. It may well be that the same mecha- 
nism that is responsible for the high- velocity absorptions at 
optical wavelengths (a density or abundance enhance ment 
in the very outer layers, see e.g. iMazzali et al.l 120 05) also 
acts to produce a NIR Mgll line that is more pronounced 
and more strongly blueshifted than usual. 



3.2 UV-optical-infrared 

While UV-optical spectra (with varying quality of the 
UV part) have been published for several SNe la (e.g . 
Kirshner et~aD Il993l : iBufano et all I20Q9I: IWang et al.l 12012 



Foley et al .1 l2012al : LMaguire et al.l l20ia ~ and combined 



optical-NIR spectra are available for many well-observed 
SNe la, complete UV-through-NIR spectra are still a nov- 
elty. In fact, to the authors' best knowledge only a sin- 
gle UV-through-NIR spectrum of a SN la (though with a 
very high quality UV part) has thus far seen the light of 
day: a max i mum- li ght spectrum of SN 201 liv presented by 
iFolev et al.l (|2Q12bh . 

The unusually complete data set of SN 2005cf in dif- 
ferent wavelength regimes ena bled us to create a set of 
five spectra combined from UV (Bu fano et alfeooil ). optical 
(jGaravini et al.ll2007h and NIR (this paper) data (Fig. [3]). 
As the spectra were typically not taken at exactly the same 
epoch, we chose the spectra closest in time and adjusted the 
fluxes of the UV and NIR spectra to match the flux level 
of the optical spectra. Our UV-through-NIR spectra which 
cover epochs from — 8 d before to — 10 d after maximum 
light, are presented in Fig. [3] 



IFolev et al.l (|2Q12bh already commented on the SED of 
SN 201 liv at maximum light being strongly peaked at opti- 
cal wavelengths. From Fig.[3]one can see that for SN 2005cf 
this is true for all epochs from 8 d before to 10 d after max- 
imum, and that the SED does not change much in this pe- 
riod. The flux blueward of 3000 A and redward of 10 000 A 
never contributes much to the total emission. The strong 
UV flux deficit in otherwise blue spectra is caused by se- 
vere line blanketing of iron-group elements at wavelengths 
shorter than - 4000 A dPinto Eastman! l2000h . The main 
route for UV photons to escape is by fluorescence into the 
optical or NIR regime, which is discussed in more detail in 
Section EL2l 
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Figure 3. Time sequence of combined UV (|Bufano et al.l BoOQ). 
optical (jGaravini et al.l [2007) and NIR (this work) spectra of 
SN 2005cf, scaled to the same peak flux. The NIR spectrum of 
the first epoch shown in the figure is an average of our —10. 2d 
and — 5.4d spectra to approximate the spectral appearance at the 
epoch of the UV-optical spectrum. 

4 MODELLING 

The discussion of spectral features in the previous section 
was solely based on a comparison of our NIR spectra of 
SN 2005cf to observed spectra of similar SNe la and line 
identifications made for those SNe. Here, we use detailed 
radiative-transfer models to discuss the spectral features in 
the NIR and the evolution of the overall SED in more detail. 

4.1 Setup 

To this end we used the Monte Carlo rad iative-transfer code 
ARTIS (|Kromer fc Siml 120091 : ISiml l2007h to obtain a self- 
consistent solution of the radiative-transfer problem and cal- 
culate synt hetic spectra for the hydrodynamical exp losion 
model W7 (iNomoto et al.lll984l ; llwamoto et "all 1 19991 ). W7 
is a one-dimensional model of a carbon deflagration in an 
accreting Chandrasekhar-mass carbon-oxygen white dwarf. 
It is a well-studied and widely-distributed model and repro- 
duces observed optical spectra and light curves of SNe la 



reasonably well (e.g . [Branch et al. 1985; Kasen et al. 2006; 
iKromer fc Simlboodfjack et al.ll20 111 : Ivan Rossum 201^). In 
the NIR, however, W7 has not been studied in great detail 
yet and the main focus ha s been laid on the ligh t curve evo- 
lution dKasen et al.l 120061 : IKromer fe Siml 120091 : Ijack et al.l 
l2Q12h . The synthetic NIR spectra of W7 presented here 
can thus be a basis for future comparisons among different 
radiative-transfer codes. 

F or our simulatio n we u sed the 'big_gf-4' atomic data 
set of IKromer & Siml ([20091 , see their table 1) with a to- 
tal of rsj 8.2 x 10 6 lines. The ejecta were mapped to a 50 3 
Cartesian grid. We then computed the propagation of 10 8 
photon packets from 2 to 120 d after the explosion. The cal- 
culation was split into 111 logarithmic time steps. To speed 
up the initial phase, a grey approximation, as discussed by 
IKromer &; Siml (|2009h , was used in optically thick cells, and 
the initial 30 time steps (i.e. the initial six days after the 
explosion) were treated in local thermodynamic equilibrium 
(LTE). Since this is well before the epoch of our first spec- 
trum of SN 2005cf it should not have any influence on the 
following discussion. 

Three snapshots of the resulting spectral evolution at 
—5.6, 1.7 and 30.9 d with respect to B-band maximum 
(which occurs 19.5 d after the explosion, at an absolute mag- 
nitude of —19.14) are shown in Fig. [4] for UV to NIR wave- 
lengths. 

4.2 Formation of the spectral energy distribution 

From our Monte Carlo simulation we do not only get the 
spectral evolution but also information on the last absorp- 
tion and emission processes. In Fig. 2] this information is 
used to colour-code below the synthetic spectrum the frac- 
tion of escaping packets in each wavelength bin which were 
last emitted by bound-bound transitions of a particular ele- 
ment (the associated atomic numbers are illustrated in the 
colour bar). White regions between the emerging spectrum 
and the colour coded part indicate wavelengths at which 
continuum processes (bound-free, free-free) contribute to the 
last emission. The coloured regions along the top of the plots 
show the distribution of photon wavelengths that bound- 
bound-emitted escaping packets had prior to their last in- 
teraction. Again the colour coding indicates the elements 
which were responsible for absorbing packets in a particular 
wavelength bin. 

Due to energy conservation the total 'area' of the emis- 
sions and absorptions, respectively, is equal. Thus the in- 
formation depicted by the colour coding in Fig. 2] can be 
used to understand how the overall SED of our model 
forms. Around maximum light (middle panel) the flux is 
strongly suppressed by line blanketing of a dense forest of 
iron-group-element lines below ~3000A. Between ^3000 A 
and ^6000 A prominent absorption features due to reso- 
nance scattering in lines of intermediate-mass elements are 
imprinted on an underlying quasi-continuum emitted by a 
dense forest of iron-group-element and intermediate-mass- 
element lines. Red ward of 6000 A the emerging flux is almost 
only due to fluorescent emission - a prominent exception is 
the Can NIR triplet. 

As visible from the inlays in the top right corners of 
the plots in Fig. [4] this behaviour is even more obvious 
in the NIR where almost no absorption occurs. The only 
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Figure 4. UV to NIR spectral evolution of the W7 model for 
three different snapshots of our radiative-transfer simulation. The 
inlays give an enlarged view to the NIR regime. Epochs are given 
relative to 5-band maximum which occurred at 19.5 d after ex- 
plosion. The colour coding indicates the elements responsible for 
both bound-bound emission and absorption of quanta in the 
Monte Carlo simulation. The region below the synthetic spec- 
trum is colour-coded to indicate the fraction of escaping quanta 
in each wavelength bin which last interacted with a particular ele- 
ment (the associated atomic numbers are illustrated in the colour 
bar). Similarly, the coloured regions along the top of the plots in- 
dicate which elements were last responsible for removing quanta 
from a particular wavelength bin (either by absorption or scatter- 
ing / fluorescence). White regions between the emerging spectrum 
and the colour-coded bound-bound emission indicate the contri- 
bution of continuum processes (bound-free, free-free) to the last 
emission. 



feature which shows some absorption contribution in this 
wavelength regime is the aforementioned feature at 10 400 A, 
which can be clearly identified as the Mgil A10926 line from 
our simulation. Although even this line has a significant con- 
tribution fro m fluorescent emission , this confirms the inter- 
pretation of I Wheeler et al.l ([ 19981 ) and justifies our use of 
the associated P-Cygni absorption trough in Section 13.1.41 
to measure ejecta velocities of Mg. Note also that contin- 
uum processes contribute to some extent to the emission in 
the NIR (indicated by the white region between the emerg- 
ing spectrum and the colour coded emission due to bound- 
bound processes). From our simulations we find that this 
continuum contribution is dominated by free- free emission. 
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Figure 5. Top panel: abundance structure (colour-coded are the 
mass fractions of elements, with the associated atomic numbers 
given in the colour bar at the bottom of the figure) and density 
profile (black line) of the W7 model. Bottom panel: comparison of 
the emissivity in the last interaction at UV (1500-3500 A; blue), 
optical (3500-9500 A; green) and NIR (9500-26 000 A; red) wave- 
lengths in the W7 model 1.7 d after maximum light as obtained 
from our Monte Carlo simulation. 

Comparing the synthetic spectrum at about maximum 
light (middle panel of Fig. U) to pre- and post-maximum 
epochs (top and bottom panels, respectively) it is clearly 
visible that the SED is dominated by different species at 
different epochs. While before maximum light intermediate- 
mass elements dominate the spectral features, they are al- 
most negligible at 30.9 d after maximum (the strong NIR 
triplet of Can is a prominent exception). Instead, the spec- 
tra are then completely dominated by fluorescent emission 
of iron- group elements. This is a consequence of the increas- 
ing dilution of the ejecta with time, which leads to lower 
optical depths allowing an observer to look deeper into the 
ejecta (which are strongly layered in the case of W7, see top 
panel of Fig. [5}. 

In the NIR this transition to an iron-group-element- 
dominated spectrum is already completed at maximum 
light, and even at pre-maximum epochs iron- group-element 
emission contributes a significant fraction to the emerging 
flux. This is further illustrated in the histograms in the bot- 
tom panel of Fig. [5] which show the regions of last emis- 
sion for different wavelength regimes at 1.7 d after maxi- 
mum light. While the NIR emissivity is strongly peaked at 
~ 5000 kms -1 , the optical emissivity has significant contri- 
butions from a wider velocity range. The UV, in contrast, 
originates from a narrow region outside the iron-group- 
element-rich inner core, since the dense forest of iron-group- 
element lines causes extremely large optical depths within 
the core so that UV photons are trapped. The substructure 
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in the histograms is a consequence of the emission by differ- 
ent species and the strong layering of the eject a. 

4.3 NIR line identification and comparison to 
SN 2005cf 

The information on the last absorption and emission pro- 
cesses presented in the last section cannot only be used to 
investigate the formation of the broad-band SED, but also 
to identify individual features in the spectra. For that pur- 
pose contributions from individual ionization states can be 
extracted from the simulation, similarly to the contributions 
of different elements shown in Fig. 2] 

We used this approach to identify individual features 
in the NIR regime of our model spectra. The results are 
presented in Fig. [6] along a time sequence of our synthetic 
spectra of the W7 model for 8 epochs between —10.2 and 
41.2 d with respect to B-band maximum. In order to im- 
prove the signal-to-noise ratio the W7 spectra at 10.2, 5.6, 
0.5 d before and 1.7 d after peak brightness were obtained 
by averaging over five, three, two and two subsequent time 
steps, respectively. This is not expected to have a signifi- 
cant influence on the physical information encoded in the 
respective spectra, as adjacent time steps are very close at 
the earlier epochs due to their logarithmic distribution. 

For comparison Fig. [6] also shows our NIR spectra of 
SN 2005cf at the corresponding epochs. In general, the agree- 
ment between our synthetic W7 spectra and SN 2005cf is 
quite good and it seems that the W7 model, designed only 
to fit optical spectra, can also represent NIR spectra at early 
epochs. That we identify many more features in the W7 
spectrum at —10. 2d than we can safely detect in SN 2005cf 
is not necessarily a flaw of W7. It could rather be due to 
the fact that these features are relatively weak (both in the 
SN and in the simulation) and we are only able to identify 
them for the W7 model because we can resolve the various 
contributions. 

Ten days after 5-band maximum the agreement be- 
tween the synthetic and the observed spectra starts to dete- 
riorate redward of 14 000 A. Although our W7 spectra repro- 
duce the characteristic Co n features which dominate this 
wavelength region at about the right epoch, their relative 
strengths do not always match those observed in SN 2005cf. 
Also there seems to be a small wavelength offset for the Co II 
15 500 A feature. This, for example, might be explained by 
uncertainties in the atomic data. Many of the lines used 
i n our simula tions which are based on the atomic data of 
(Kurucz 2006) are only theoretically predicted so that both 
oscillator strengths and resonance wavelengths could be at 
fault |j or so me w eaker lines may be m issing completely. 
iKasenf (|2QQ6h and iKromer fc Siml ([20091 ) have shown that 
this can have a significant influence on NIR observables. 
Moreov er, the simplified exc itation treatment used by AR- 
tis fsee lKromer &; Simll2009l for a detailed description) may 
lead to inaccurate populations of excited levels, particularly 

6 For example, for the Coll lines given in Table 0] we com- 
pared the wavelengths found in the atomic data base of IKuruczl 
( 2006) with the corresponding level energy differences in the NIST 
database (http://www.nist.gov/pml/data/asd.cfm). From this we 
find typical offsets on the order of <300kms _1 . 



Table 4. Properties of 8 Coll lines that dominate the X-band 
spectrum of SNe l a a mon t h afte r maximum light, based on the 
atomic data set of IKuruczl (|2006l V 



Wavelength 


Upper level 




Lower level 




(A) 


Conf. 


Term 


J 


Conf. 


Term 


J 


20 913 


3d 7 ( 4 F)4p 


z 5 F 


1 


3d 6 4s 2 


a 5 D 


1 


21211 


3d 7 ( 4 F)4p 


z 5 F 


2 


3d 6 4s 2 


a 5 D 


2 


21351 


3d 7 ( 4 F)4p 


z 5 F 


4 


3d 6 4s 2 


a 5 D 


4 


21509 


3d 7 ( 4 F)4p 


z 5 F 


3 


3d 6 4s 2 


a 5 D 


3 


22 209 


3d 7 ( 4 F)4p 


z 5 F 


5 


3d 6 4s 2 


a 5 D 


4 


22 482 


3d 7 ( 4 F)4p 


z 5 F 


2 


3d 6 4s 2 


a 5 D 


1 


23 619 


3d 7 ( 4 F)4p 


z 5 F 


3 


3d 6 4s 2 


a 5 D 


2 


24 603 


3d 7 ( 4 F)4p 


z 5 F 


4 


3d 6 4s 2 


a 5 D 


3 



at late epochs when non-LTE effects become increasingly 
important. 

An interesting result, which will be unaffected by the 
aforementioned uncertainties on a qualitative level, is the 
fact that starting two weeks after maximum the emission in 
the various NIR bands is dominated by different iron-group 
species. While the J-band is dominated by Fen emission 
the H- and if-band emission originates mainly from Coll 
(see Figs. [4] and [6}. Since the ionisation potential of both 
species is of the same order, the J — H and J — K colours 
are potential probes to the Co/Fe-ratio in SNe la ejecta. 
Compared to the observed spectra of SN 2005cf, the W7 
spectra have too much flux in the H and K bands relative 
to the J band at these epochs (Fig. [6]). This might indicate a 
deficiency of stable iron- group elements in the ejecta of W7. 
However, on a quantitative level the uncertainties in the 
atomic data and the radiative-transfer treatment currently 
prevent any firm conclusions on this issue. 



4.4 Constraining the extension of the iron core 

Like most other features in the NIR spectra of SNe la the 
four characteristic Coil humps that dominate the if-band 
spectrum a month after maximum light are blends. How- 
ever, our simulation shows that in this case the number 
of involved transitions is comparatively small. In fact, the 
main contribution comes from no more than 8 Coll emis- 
sion lines, which are characterised by particularly low-lying 
upper levels. Apparently, the population of other Coll lev- 
els from which if -band transitions can originate is very low. 
The wavelengths and configurations for the 8 Co II lines are 
given in Tabled 

We note that the small number of contributing lines 
makes it possible to deblend the if-band region and derive 
the extension of the emitting Co-rich core independently 
from invoking a complex modelling machinery. As a simple 
recipe, this can be accomplished by fitting Gaussians to each 
contributing emission line, keeping the central wavelengths 
fixed and enforcing a common full width at half maximum 
(FWHM) in velocity space. The FWHM which - besides the 
individual intensities - is a fit parameter, then gives an es- 
timate for the extension of the emitting Co-rich core. Note, 
however, that there is not necessarily a one-to-one corre- 
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Figure 6. NIR time sequence of our synthetic spectra of the W7 model, in comparison with the corresponding NIR spectra of SN 2005cf. 
An identification of prominent lines in the W7 spectrum as obtained from our Monte Carlo simulation is also shown. 



spondence between the two quantities so these velocities are 
only indicative. 

Applying this technique to our spectra of SN 2005cf us- 
ing 8 Gaussians corresponding to the aforementioned Coll 
lines (Table |4|, we find that the wavelength region be- 
tween rsj 20 800 A and ~ 25 000 A is adequately reproduced 
with a FWHM velocity of ~ 11 000 km s -1 for the spectra 
at +31 and +41 d. Slightly smaller FWHM velocities of 
~ 10 000 to 10 200 km s" 1 are derived for the +29 d spectrum 
of SN 2002bo and the +32 and +42 d spectra of SN 1999ee. 

These numbers are in rough agreement with the ex- 
tension of the iron-group-element zone in W7 where those 
elements are abundant at the 10 per-cent level at a velocity 
of 10 600kms _1 , and drop to a mass fraction of 1 per cent 



at 12 300kms _1 . The fair match could already be expected 
from the rough agreement between the spectra of SN 2005cf 
and the model spectra of W7. It is, however, no validation 
of the W7 model, since in detail there are numerous other 
differences between the model and the observed spectrum. 



5 CONCLUSIONS 

We have presented NIR spectroscopy of SN 2005cf covering 
epochs from — 10 d to +42 d with respect to B-Band maxi- 
mum. Together with archival data at other wavelengths this 
makes SN 2005cf one of the currently best-observed SNe la 
and allowed us to compile combined UV-optical-NIR spectra 
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at five epochs between ~ 8 d before and ~ 10 d after maxi- 
mum light. In this respect SN 2005cf is a unique object since 
to our best knowledge thus far only a single UV-through- 
NIR sp ectrum of a SN la has been published by I Foley et al.l 
(|2Q12bh for SN 2011iv. 

A com parison of our spe ctra to observations of 



SNe 1999ee (lHamuv et al.l l2QQ2allbh , 2 02bo (jBenetti et all 
l2004h and 2QQ3du (jStanishev et all l2007h shows that 
SN 2005cf is a perfectly normal SN la at NIR wavelengths. A 
similar con clusion has already be en drawn at optica l wave - 
leng ths bv iGarayini et al.l (|2QQ7h . IPastorello et alJ (I2007T ) 
and lWang et al.l (l2009l ~ 

We have also performed radiative-transfer simulations 
to obtain NIR spec tra for the standard hydrodynamical ex- 
plosion model W7 ([Nomoto et al.|[l984h . Comparing the ob- 
tained synthetic spectra to our spectra of SN 2005cf we find 
that W7 reproduces the key properties of normal SNe la 
not only at UV-qptical as reported in earlier work (e.g . 
Branch et a~lll985l : iKasen et al1l2006l : iKromer fc Sim| [2QQ9: 



Foley et al.l 



2012b|) but also at NIR wavelengths. From a 



detailed analysis of the emission processes in the radiative- 
transfer simulation we have identified the main spectral fea- 
tures in the NIR. We find only a single clear P-Cygni line: 
the Mgll A10926 triplet which is present at early epochs. 
Apart from that the NIR part of the spectrum is almost 
devoid of absorption features. Instead, fluorescence is the 
driving mechanism that shapes the SED of SNe la. Before 
maximum light both intermediate-mass elements and iron- 
group elements contribute to the fluorescent emission. After 
maximum light almost all the flux originates from iron-group 
elements, though different iron- group species contribute to 
the various NIR bands. This makes the J — H and J — K 
colours a few weeks after maximum light potential probes 
to the ratio of different iron- group species. While most of 
the fluorescent emission is strongly blended, we find that 
the characteristic Co n humps which form at about a month 
after maximum light in the K band result from a series of 
relatively unblended lines. This provides a simple possibil- 
ity to constrain the iron-group-element-rich regions of the 
eject a from observed spectra. 
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